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1 Introduction and Background Material 
1.1 Introduction and Summary 

In recent years the concept of squeezed state has become central in quantum 
optics both from the theoretical and experimental point of view In the 
simplest one-mode case a squeezed state is defined here as a displaced squeezed 
vacuum of the form |a,C) = D{a)S{()\0) obtained by applying the squeez- 
ing operator S{Q = e^'' " -'/^ a and the displacement operator D{a) = 
gQa^-a a ^^le vacuum state |0) (i.e. the ground state of the one-dimensional 
harmonic oscillator). ( ~ re*^^ and a — \a\e'^^ are complex numbers, and 
a — {q + ip) / V2, = (9 ~ ip) I in terms of the normalized coordinate and 
momentum operators q and p, where = il. The term "squeezing" comes 
from the fact that for cos29 > tanhr(> 0) the dispersion Aq of the coordinate 
operator in the squeezed state is smaller than the vacuum state value 1/^/2 : 
the uncertainty of the value of q is squeezed compared to the vacuum value. 
If one performs a rotation in the gp-planc ( "phase space" ) to new canonically 
conjugate operators qe — cos 9 q + sin 9p, pg ~ ~ sin 9 q + cos 9 p the dispersions 



*This is a reproduction, with minor corrections and an added Appendix, of a paper pub- 
lished in J. Bertrand et al.(Eds.), Modem Group Theoretical Methods in Physics, Kluwer 
Academic Publishers 1995. 

^We follow the convention used e.g. in the first and third references in sometimes the 
squeezing operator is defined with the opposite signs in the exponent. 
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Age and Ape will be equal to e^^ /^/2 and e''/-\/2, respectively, and the covari- 
ance /S.[qepe) □ is zero. This implies that the variance matrix of the original set 
qp {a. symmetric 2x2 matrix with the variances (Ag)^ and (Ap)^ as diagonal 
elements and A(gp) as off diagonal elements) has its determinant equal to the 
minimum allowed value j. This can actually be taken as a characteristic prop- 
erty of squeezed states: they are the states (amongst a priori even mixed states) 
that give equality in the relation (Ag)^ (Ap)^ — [A(gp)]^ > j (sometimes called 
the Schrodinger-Robertson uncertainty relation). The coherent states {( — 0) 
are characterized by having the variance matrix diagonal with diagonal elements 
equal to ^ . An equivalent characterization of squeezed states uses the fact that 
I a, C) is the (up to a phase factor unique) normalized vector annihilated by the 
transformed annihilation operator b — coshr {a — a) + e*^^ sinhr (a^ — a*). 

The definition of squeezed states we prefer is as the coherent states (with 
vacuum as isotropy vector) of a representation (rep) of a group Hx M(2,R), 
which is a semi-direct product of the 3-dimensional Heisenberg-Weyl group H 
and the 3-dimensional metaplectic group M(2,R), a two-fold covering of the 
group SU(1,1). 

In experimental realizations of squeezed states the state often corresponds 
to a two-mode squeezed state generated by applying the two-mode squeezing 
operator 82(0 — a2-(a-^ a^) ^j^g vacuum state. 

The purpose of this paper is to derive the spectrum and the (generalized) 
eigenfunctions, in three different forms, of the squeezing generators (a^^ — a^)/4i 
and (a| — ai a2)/2i, and also of another quadratic (in a and a^) generator 
introduced in the next section. Using these results the spectrum and eigenfunc- 
tions of the most general second degree polynomial hermitian Hamiltonian (not 
necessarily bounded below) can be obtained. 

The problem of deriving the eigenfunctions of the squeezing generators has 
been studied before ||], but the results were not complete, in particular the 
doubling of the spectrum in the one-mode case was not obtained. 

In the next section the group theoretical background of the squeezed states 
will be quickly presented, and after that the three different state representa- 
tions will be given: the photon number (harmonic oscillator excitation number) 
representation (the n-representation) . the ordinary configuration space repre- 
sentation (the q-representation) and the Fock-Bargmann holomorphic function 
representation (the z-representation) . The last representation has the advantage 
that the generalized eigenfunctions will be holomorphic functions (but of course 
not normalizable) , whereas in the ^-representation we expect to get tempered 
distributions, and in the n-representation tempered sequences (growth bounded 
by some power of n). In general the starting point will be the rt-representation, 
where the method of an earlier paper ||] can be used. 

2We use the notation A{AB) = {{A' B' + B'A')/2), where A' = A - (A) etc; A{AA) = 
(AA)2. 
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1.2 Squeezed States as Coherent States of HxM(2,R) 

The presentation given here is essentiaUy weh-known Q , so we only sketch the 
derivations. The condition for a hnear inhomogeneous transformation of {a, a'} 

b — iia + i'a'+a, = v* a + fi* a' + a* (1) 

to give a new set {b, b^} satisfying the canonical commution relation (CCR) 
is l/xp — — 1, which means that the homogeneous part of the transforma- 
tion is an element of the group SU(1,1). It is known that any two irreducible 
representations {irreps) of the CCR are unitarily equivalent. Now the original 
rep by {a, a^} is clearly irreducible, and it is not hard to convince oneself that 
this also holds for the rep by {b, 6^} — one can show that there is a unique 
"new vacuum" |0)b with 6|0)b = 0. This means that one gets a set of unitary 
operators {V{a, g) : a £ C, g £ SU(1, 1)} , determined up to phase factors, 
such that 

b^V{a,g)^aVia,g) (2) 

Composing two successive transformations and using irreducibility one realizes 
that one has a rep up to a phase factor of the group ISU(1,1), "inhomogeneous 
SU(1,1)", with the composition {a\g') ■ {a,g) — {a' + g'a,g'g), where g'a = 
n'a + v' a* . It is then known that the obstruction against choosing the phase 
factors of V{a, g) in such a way that one gets a true rep is present already in 
the 2-dimensional translation subgroup {a}: since a and don't commute one 
has to make a (central) extension of this group to the 3-dimensional Heisenberg- 
Weyl group H= {{a, s) : a £ C, s G R} with the composition (a',s') (a, s) = 
(a' + a, s' + s + Im[a'a*]). Furthermore it turns our that the resulting rep of 
SU(1,1) is two- valued, so that if one introduces the two-fold covering of SU(1,1), 
called the metaplectic group M(2,R), one obtains a true unitary rep of the group 
Hx M(2,R). Here already the restriction to the group H is irreducible, whereas 
the restriction to M(2,R) gives the sum of two irreps. Corresponding to a 
factorization of a general group element in the form 

(a,5,s) = (a,e,0)-(0,g,0)-(0,e,s) (3) 

we can write a general representing unitary operator in the product form 

U{a,g,s)^D{a)U{g)e'\ (4) 

The standard Perelomov definition of a set of coherent states of a given rep of 
a group G starts by choosing a particular isotropy vector of the rep which is 
invariant (up to a phase factor, i.e. gives a one-dimensional rep) under a suitably 
chosen subgroup H of the given group. The set of corresponding coherent states 
is the set of states (ignoring differences in phase factors) obtained by applying 
all unitary operators in the rep on this vector. The set can then be parametrized 
by the elements of the homogeneous space G/H. 
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In our case we choose the vacuum state as isotropy vector (as a more general 
case one could choose an arbitrary number state |n) to get displaced squeezed 
number states), and we then want to factorize the SU(l,l)-representing operator 
as 

f/(ff) = 5(C)-e*2^-^° (5) 

where as earlier 5(C) = e(^*«'-^"*')/2, and Jo = {aU + aa^)/4 = {N + i)/2 is 
the generator of a U(l) subgroup of SU(1,1). As 

5(0^0 5(0 =coshr- a- e'^^sinhr- at, e-'^*'"^" a e"*'"^" = e'^ a, (6) 

we get by comparing with U{g)^a U{g) = + ua^ the relations 

^ = e*"^ cosh r, -e^^^^-v) ginh r. (7) 

Jo together with Ji = {a^"^ + 0^)/^ and J2 = (a'''^ — a^)/4i generates a rep 

[Jo,Ji]=iJ2, [Jo,J2] = -iJi, [Ji,J2] = -iJo (8) 

of the Lie algebra su(l,l) with Casimir operator C2 = Jq — Ji — Jf = 
We get two irreps of the two-fold covering of SU(1,1) belonging to the discrete 
series with k = j and | (C2 = k{k — 1) in an irrep), corresponding to 
the spectrum SpJo = (spanned by the even number states), and 

{f 5 i> • • •} (spanned by the odd number states). Wc shall also study the gener- 
ator = Jq + Ji = (a + at)^/4. Jq, J2, and generate reps of the three 
different types of one-parameter subgroups of SU(1,1): elliptic, hyperbolic, and 
parabolic, respectively. Since two irreps are involved, from general results it 
follows that the spectrum of J2 is twice the whole real line, and for K+ twice 
the positive half-axis. It remains to give explicit exprc!ssions for the correspond- 
ing generalized eigenfunctions. As mentioned earlier this will be done in three 
different forms, which we now proceed to define. 

1.3 Three Different Representations of Oscillator States 

The connection between the three different representations introduced earlier, 
the n-,q- and ^-representations, can be made via the ON bases correspondences 
as follows: 

{\n)} ^ {un{q) = NnHniq) ■ e-«'/2} ^ where ^„ = (0F2"n!)-i 

(9) 

\'^) =Yl(^ri\n) ^ ifi{q) ='^anUn{q) ^ (j){z) = ^anz'' /\/n\ = el^^l '"^ {z*\lp), 

(10) 

where {a\'f) is the so-called coherent state representation; 

ll<^ir = E = / l'^(9)r dq =n-' j J \ct>{z)\'e-\^\' dH. (11) 
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The generators of the representation of our group has the forms 

a ^ {q + d/dq)/V^ ^ d/dz, ^ (q — d/dq)/y/2 ^ z; 
Jo : (flta + aat)/4 ^ {-d'^/dq^ + q^)/A = Hoscf^huj ^ {zd/dz + l/2)/2; 
Ji : {a^^ + a^)/4^ {dydq^ + q^)/A^ {d^dz^ + z^)/4; (12) 

J2 : (at2 - a^)/Ai ^ i{qd/dq + l/2)/2 ^ i{d^/dz^ - z^)/4; 

K+ : (at + a)V4 <-> gV2 ^ (dVdz^ ^ 2zd/dz - z^ + l)/4. 

We shall need formulas taking us from (p{q) to 0(2;) and vice versa (Bargmann 
and inverse Bargmann transforms). These are easily obtained using the integrals 
for a„: 



u^{q)ip{q)dq^TT-' j J {z"' /Vn\)4>{z) e-^'^ d' z (13) 

and the generating function for the Hermite polynomials 

exp{2yt - t') = ^ H,,{y) e/n\. (14) 

The result is 

cPiz)^ J K{z,qMq)dq; ^{q) ^ j j K{z\q)<i>{z) e'^'^" d^ z, (15) 
where the kernel 

K{z, q) = TT-i/'* exp(-zV2 - + V^zq). (16) 

2 Generalized Eigenfunctions of the Squeezing 
Generator J2 

We start by deriving the eigenfunctions in the ^-representation. Since the eigen- 
functions belong to a continuous spectrum, the sequence {a„} we obtain will not 
be square summable, only power bounded (i.e. a tempered sequence, dual to 
rapidly decreasing test sequences). Se e.g. j^] for details about the mathematical 
background. So we have to solve the eigenvalue equation 

J2IA) = [(at2-a2)/4*]|A) = A|A). (17) 
With |A) = ^ fn{^)\n) this leads to the recursion formula 



v/n(n-l)/„_2(A) - v/(" + l)(ri + 2)/„+2(A) = ^4A/„(A), (18) 

evidently connecting only the even or the odd number states, corresponding to 
a doubling of the spectrum. To get rid of the square root we put {N[\) is a 
normalization constant) 



/„(A) = N{\)^T{[n + l]/2)/r(n/2 -f 1) a„ (19) 
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to get the relation 

(n + l)a„+2 + j4Aa„ - na„_2 = 0. (20) 

We first treat the case of even numbers n and put n = 2m and bm = a2m, 
so that ^ 

{m + -)brn+i + i2Xbm - mbrn-i = 0. (21) 

This recursion relation can be solved by the standard Laplace method intro- 
ducing the expression b^. — J t"^ y{t) dt in the relation, which gives after partial 
integrations 

j t"'[-d{fy)/dt + ty/2 + i2Xy + dy/dt] dt - If"(l - t^)yit) = (22) 

where I... gives the contribution from the ends of the integration path. Taking 
this path to be from -1 to 1 (so that the corresponding contribution vanishes) 
and solving the simple first order differential equation for y{t) we get after an 
integration variable change leading to an integral from to 1, and using the 
standard integral representation for the hypergeometric function 2^1 (which for 
convenience we denote by F) the result 

-.„. = (-l)"2-V^!Wl^±ffilZi±I^f 1/4 - .A, 1/2; 2) ,23) 

One could also have found this result directly by comparing the recursion 
relation with a suitable contiguity relation for the hypergeometric functions Q . 
A transformation formula for transforming z — > z/(z — 1) in gives 

F(-m, 1/4 - a, 1/2; 2) = (-1)"F(-to, 1/4 + iX, 1/2; 2). (24) 

Together with the relation, valid for real A, 

F(-m, 1/4 - iX, 1/2; 2)* = F(-m, 1/4 + iX, 1/2; 2) (25) 

this implies, using the series expansion of F , that F is for even m a real, 
even polynomial of degree m in A, whereas for odd m it is i times a real, odd 
polynomial of degree m. With a change of definition of the normalization factor 
we can write 

/2m(A) = Ar(A) V(l/2)™/™! F{~m, 1/4 + iX, 1/2; 2), m = 0, 1, .. (26) 

using the notation {a)m = r(a + rn)/T{a). 
Treating the odd n case in a similar way we obtain 

/2™+i(A) = A^'(A)V(3/2),„/m!F(-m, 3/4 + tX, 3/2; 2), m = 0, 1, ... (27) 
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We now have to prove the (generahzed) orthonormahty and completeness 
relations of these generalized eigenfunctions, and in particular determine the 
normalization constants. Although this can be done using (26) and (27) and 
suitable integral representations it is simpler to observe that we have obtained 
special cases of the Pollaczek polynomials Q u 

Pm{\, b) = i"'^i2b)m/m\F{-m, b + ^A, 26; 2), m = 0, 1, .. (28) 

which for every 6 > form a conmlete orthonormal set of polynomials on the 
real line with the weight function □ 

Pb(A) = 22^-i|r(6 + zA)|V^r(26). (29) 

If we choose positive normalization factors 

^(A)-y^^M^, N'{X) = y/p3/4(A) (30) 

and define 

|A, e)=J2 /2m(A)|2m), /^^ = (-*)™^Pi/4(A) P™(A, 1/4), (31) 

|A,o) -^/2™+i(A)|2m + l>, h,n+i - (-*)'"^P3/4(A)P™(A,3/4), (32) 
then the generalized orthonormahty relations would be 

(A,e|A',e) (A,o|A',o) =(5(A- A'), (A, e|A', o) = 0. (33) 

Whereas the last relation is obvious, the first two just express the completeness 
of the corresponding set of Pollaczek polynomials, and thus hold true. 

In a similar way the orthonormahty conditions for the Pollaczek polynomials 
imply the completeness relation for the generalized eigenfunctions, since this 
relation says that 

J[\{X,em' + \{X,om']dX (34) 
for an arbitrary state |'(/') = ^ c„|n), which can also be written in the form 

'^\Cn\'^ = ^[clmC2n / /2m ( A)* /2n ( A) dA + C2„i+i C2„+l / /2m+l (A)* /2«+l (A) rfA] . 
m,n 

(35) 

^This definition difi'ers with the factor from the definition in the original article, and 
makes Pm a real polynomial with positive coefficient for the highest power term 
■'see Appendix 
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As Po{X,b) = 1 it follows from (31) and (32) that 



(0|A, e) = y'pi/4(A), (1| A, o) - ^P3/4(A) (36) 

To derive the form of the eigenfunctions in the z-representation one could of 
course start from the corresponding eigenvalue equation 

i-dydz^ + z^)^xiz) = z4A0a(z). (37) 

We shall instead use the result we have obtained in the n-representation and 
apply a standard integral representation for the F function to obtain from (26) 



Jo 



3/4-jA^-[^_^-j-3/4+jA 



= A(A)e-^'/2 / i-3/4-'^(l-f)-3/4+'Ag.^trf^^ where^(A) = y^Pi/4(A)-— 



r(l/2) 



Xl/4 + zA)|2- 
(38) 

Use of a standard integral representation for the confluent hypergeometric func- 
tion iFi results in the final form 



cl>x,e{z) = ^pi/4(A)e-^'/2 ^F,{l/4: ~ i\ 1/2; z^). (39) 
In a similar way we obtain 

<l>xA^) = 7p3/4(A) e-^'/2 z iF^ii/A - i\ 3/2; z"). (40) 

As the second parameter of the iFi function is 1/2 or 3/2, respectively, the 
functions can be expressed in term of the parabolic cylinder functions Du{\/2z) 
and Di,{-V2z) with i/ = -1/2 + i2\ (see §], p.ll7). 

Finally we shall derive the form of the eigenfunctions in the g-representation. 
In this case we shall solve the corresponding eigenvalue equation, which can be 
written 

qd^x{q)ldq= -{- + i2\)^x{q). (41) 

We have a singularity at g = and can choose independent solutions on the 
positive and negative half-axes. These solutions can be put together to give 
solutions even and odd, respectively, in q, corresponding to even and odd number 
solutions. We then get 

^UD = (2-)-^/^^T7W' = (42) 

The factors (27r)~^/2 are introduced to ensure the generalized orthonormality 
relation 

^y{qrMq)dq={27T)-'2 / |g|-l+*2(V-A)^g 
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= (27r)-W exp[i{X' - X)t]dt = S{X' - X). (43) 

J —oo 

It is clear that ip\e{q) and ipxo{q) should correspond to the earlier determined 
eigenfunctions up to A-dcpendent phase factors. To determine these we calcu- 
late the scalar product with the vacuum and 1-photon state, respectively, and 
compare with (36): 

7r-i/4 J e-'^'^'ipxeiq) dq = 7r-i/4(2/7r)i/2 ^-i/2-»2Ag-^V2 ^ (27t')-^/^^ 

(44) 

--^/^V2 J qe^^'/\M dq = 2(2.3)-l/4 r(3^pA) ^ ^^^^ 
From this we finally get 

(g|A,e) =exp[z{argr(l/4 + a) + Aln2}](27r)-i/2^^i^, (46) 
{q\X, o) = exp[z{argr(3/4 + »A) + A In 2}] (27r)-^/^ |g|i/^4!2A " (4^) 

It should be observed that although the eigenfunctions in the ^-representation 
are ordinary functions (though with a singularity at the origin) they are not 
square integrable but should be considered as tempered distributions on the real 
line, similar to the eigenfunctions e*^'^ of the momentum operator p — —i d/dq. 
See e.g. Q for more details on the mathematical nature of the generalized 
eigenfunctions of the SU(1,1) generators. 



3 Generalized Eigenfunctions of the Generator 




We have seen in (12) that the parabolic generator has the forms (a^ -|-a)^/4, 
g^/2, and [d"^ /dz^ + 2zd/dz — z'^ + l)/4 in the different representations. As in 
the previous case one could of course start from the eigenvalue equation in the 
n-representation and solve the corresponding difference equation, which again 
involves either even or odd number states, corresponding to a doubling of the 
spectrum. In this case it is simpler to use the form of the operator in the 
g-representation, where the eigenvalue equation 

(gV2)<^r,(g) = V Vniq) or (<z' - 2ry)^^(<z) = (48) 

shows that the spectrum of is the half axis [0, oo) and that we can choose 
as properly normalized even and odd solutions 

ip,,M = {ll2^)[5{q-^)+5{q+^)], (49) 
ip^o{q) = {l/2^)[5{q-^)-5[q+^% (50) 
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where the normahzation is determined (up to a 77-dependent phase factor) by 
demanding 

j Vr,'e{q)*V'ne{q) dq = S{r]' - 77) etc. (51) 

Completeness of the set of eigenfunctions follows from the fact that together the 
two families give a family of 5- functions covering the whole g'-axis. 

The corresponding eigenfunctions in the n- and ^-representations are ob- 
tained from (13) and (15) as 

\r],e) = '^g2m{v) |2m), |ry,o) = ^g2m+iiv) |2to+ 1), where (52) 

g2m{v) = j Vr,e{q) U2m{q) dq = N2m{'i-/V^) e~'^ H2m{\/^) (53) 

g2m+l{v) = j V-qo{q)u2m+l{q)dq = N2m+l{^l\/%)e~'^H2m+l{\J%) (54) 

^e(^) = j K{z, q)vM dq = 77-^4(1/^) g-^e-^'/^ cosh(2^^), (55) 

f^o{z) = j K{z, q)vM dq = t,-^I\\I^) e-'^e-'"''^ sinh(2V^z). (56) 

4 Spectrum of the General Second Degree Hamil- 
tonian 

A general hermitian second degree polynomial in {a, a)} (or in {q,p]) can be 
written 

P2(a, a)) = A{a)a + a a^) + Be^^a^ + Be-'^'a)^ + Ce**a + Ce-**a+ + DI 

= 4[AJo + B{cos $ Ji + sin $ J2)] + C%/2(cos * g - sin 5' p) + D /. (57) 

Since we have a real linear combination of a set of generators of the group 
HxM(2,R) we can conclude that if we have a Hamiltonian of type (57) with 
arbitrary time-dependent coefficients A{t), B{t), C{t) and D{t), the correspond- 
ing time development will be given by a unitary operator of the form (4) , and 
hence a squeezed state will always (up to a phase factor) develop into a squeezed 
state. This clearly also holds for the set of displaced squeezed number states for 
a given n. 

We shall now determine the various possible spectra for the operator (57). 
It is no essential restriction to assume A, B, and C non-negative, and D = 0. 

HA ^ B we can transform away the linear terms (in C) and get the 
quadratic terms to a form proportional to Jo (if A > B) ot to J2 (if A < B) 
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by a substitution a ^ iJ,a + + a, which according to (2) can be achieved by 
a unitary transformation. We get the results 

UA> B: 4VA^ -B'^Jo- [C'^/2{A'^ - B^)][Acos{<^> - *) - Bcos^-J^SS) 
UA<B: - J2 - [C'^/2{A'^ - B'^)][Acos{^ - - B cos*] (59) 

In both cases we should choose a = -C [Ae"'* - Be'(*-*)]/[2(A2 - B^)]. 
If A = B the phase rotation a — > e~'*/^a results in the form 

2Aq'^ + C\Pi cos(* - $/2) g - C\/2 sin(* - $/2) p. (60) 

Thus if C = 0, or C ^ and $ = 2^- (mod 27r) we get the form 

2yl(g±C/2^/2^)2-CV4A. (61) 

Otherwise, i.e. if C 7^ and $ ^ 24" (mod 27r) a translation in q leads to the 
form 

2A{q^ - (C/V2 A) sin(4' - $/2) p] - {C'^j^A) cos2(* - $/2). (62) 
From the relation 

exp(igrV36)(9^ + bp)exp{-iq^/3b) = bp (63) 

it follows that in this case the operator can be transformed to the form 

-C\/2sin(*-$/2)p (or x q). (64) 

Prom (58), (59), (61), and (64) we can read off the different possible spectra 
of p2(a,a^) in (57). A = B is evidently a critical point. Here we can have as 
spectrum either a half-axis taken twice or the whole real axis. This situation is 
not stable: a small perturbation can lead either to a discrete equidistant semi- 
bounded spectrum (if A > B) or the whole real line taken twice (if A < B). 

5 The Two-Mode Case 

In the two-mode case (two-dimensional harmonic oscillator) we get an infinite 
set of irreducible (true) representations of SU(1,1) by putting 

Jo = {a\ai + a\a2 + l)/2, J\ = {a\al + «i«2)/2, J2 = {a\al - aia2)/2i. (65) 

For the Casimir operator we obtain 

C2 = Jo^- Jf- j| = ([A7V]2-l)/4, AN = Ni - N2 = alai - ala2. (66) 

Since the spectrum of AN is the set of all integers, we get for the spectrum of 
C2 = k{k — 1) as possible k values 1/2, 1, 3/2,.., i.e. all the irreducible SU(1,1) 
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reps in the discrete class . All except the first one (with An=0) occur twice. 
It is evidently enough to study those with An > 0. We get as spectrum for 
Jo {(An + l)/2 + n2; n2 = 0, 1, 2, ..}. The generalized eigcnfunctions for J2 and 
if+ can be derived in the n-representation (basis {|ni,n2)}) and from there 
obtained in the z- and g'-representations (functions (p{zi,Z2) and (fi{qi,q2)). We 
shall only give the results for both operators in the first two representations. The 
integrals leading to the (/-representation seem rather intractable in the general 
case; we only give the result for J2 in the case An = 0. 

The results for the two- mode squeezing generator J2 (An = 0, 1, ..) are 

00 

|A,An) = ^/„(A,An)|n + An,n), UX, An) = ^/^) Pn{X,c) (67) 


.^A(^i,^2;An) = VPc(A)/r(2c) " exp{-ziZ2) iFi(c - iX, 2c; 221^2). (68) 

P„(A,c) and pc{X) are defined in (28) and (29). Here c = (An + l)/2. 

For An < one should substitute |An| for An in (67) and (68), interchange Zi 
and Z2, and change |n + An, n) to |n, n + | An|). 

The generalized orthonormality and completeness properties can be written 

(A', An'IA, An) = Sau'Au HX' - A), (69) 
/ |A, An)(A, An| dX = /, identity operator. (70) 

An 

In the ^-representation we only give the result for the case An = 0: 



</?A(gi,g2;0) 



1 91-92 



7r|r(l/2-HiA)| 



qi + 92 



Kix{\ql-ql\/2), (71) 



where K^, is the modified Bessel function of the third kind. 
For the two-mode generator = Jq -\- we get 



|7?,An) =^5n(»?,An)|n+An,n), gn = (-l)"%/2e-''(27?)^"/y L^^ 

(72) 

where L" is the generalized Laguerre polynomial. 

(j>rj{zi, Z2; An) = \/2e-''(0i/02)^"/' e-'^'^'lAniV^V^), (73) 

where I,y is the modified Bessel function of the first kind. 

Changes for An < 0, orthonormality and completeness relations are analogous 

to those in the previous case. 
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6 Dedication 



This paper is a contribution to professor Guy Rideau on the occasion of his 
retirement. Rideau has been an important member of the group of scientists in 
Paris and Dijon that over the years has inspired my interest in group represen- 
tation theory. May his retirement be as pleasant as mine! 

A Completeness and orthonormality of the set 
of polynomials {Pn(A, 6); n = 0, 1, 2..} 

In this Appendix we want to show that the set {P„(A, b);n — 0,1, 2..}, where P„(A, b) = 
i"-\/(26)„/n! 2-Fi(— TO, b + iX, 2b; 2), satisfies the relations 

/oo 
P,„(A, b) Pn (A, b) pb{X) dX = Smn, (orthonormality) (74) 
-OO 

OO 

P„(A, b) Pn{X',b) = S{X - X')/pb{X), (completeness) (75) 

n=0 

where pb{X) = 2^''-^\r{b + iX)\^ /7rT{2b). 

These two relations can be proved using appropriate integral representations 
for the 2P1 function; we shall do this, but first make an alternative derivation 
of the completeness relation, based on the recursion relation satisfied by the P„ 
polynomials. This derivation has the advantage of giving the expression for the 
weight function p. 

Using the definition of P„ and the contiguity relation for 2 Pi referred to in 
[5] we find the recursion relation generalizing (18) 



c„_iP„_i + c„P„+i = A P„; c„ = -^{n+l){n + 2b); Pq = 1. (76) 
Using the notation P^ = Pn{X' , b) we can derive from ( [76| ) 

M 

CM^M = {X-X')Y,PnPn; Am = P'mPm+i - P'm+iPm ■ (77) 

ri=0 

This is the difference equation analogue of a well-known formula for linear second 
order differential equations. We shall derive the leading asymptotic behaviour 
of Am as M — > 00. To do this we need the corresponding asymptotic behaviour 
of the 2 Pi function. 
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A.l Asymptotic behaviour of 2-Fi(— ri, 6 + zA, 26; 2) as tt, — > cxd 

Using a standard integral representation of the 2F1 function and dividing the 
integration interval into parts (0, 1/2) and (1/2, 1) we can write, with N{X, b) = 
T{2h)/T{h + i\)T{b^i\), 

2Fi(-n,6 + zA,26;2) = 7V(A,6) / t^+'^-\l ~ - 2tYUt 

Jo 

= N{X, b) f^'\t{l - t)]"-! I (^^y + (-1)" (^^) "''I e"ln(l-2t) 

(78) 

It is clear from the form of the last integral that the dominant contribution 
for large n comes from the behaviour near 0. Introducing x — — nln(l — 2t) 
as new integration variable and expanding to the first two leading orders in 
x/n (actually the full asymptotic expansion can in principle be obtained by 
continuing this expansion), we get after some calculations, using among other 
things the standard integral representation of the gamma function, the result 
(we assume A real) 

2Fi(-n,6 + iA,2&;2) = (-i)"[2r(26)/|r(fe + a)|](2n)-''{cos(Aln2n - ^ - n7r/2) 

- (1/n) [bV\FT>?] cos(A In 2n~ipi- rm/2) + 0(l/n2)f79) 

where ip = argr(6 + iA), ipi = argr(6 + l + ~ (f + arctan(A/6). 

( [79I) is the result we shall use to derive the asymptotic behaviour of Aa/ . 



A. 2 Asymptotic behaviour of A^/ as M — > 00 and the com- 
pleteness relation for the set {Pn} 

Using the result (|79| ) in the expression Am — P'j^^Pm+i — P'm+iPm we obtain 
after some elementary calculations, using also the relation T{M + b) /T{M + 1) w 

Am = [r(2&)/22''-2|r(6 + iA)p](l/A/){sin[(A - A')ln2M] + (A - A')0(1/A/)}. 

(80) 

From ( pO| ) we finally obtain the completeness relation of the polynomials 
{Pn} in the form 



V P„(A, b) P„(A', b) = lim cmAm/(A-A') = [',iT{2b)/2"'-^\T{b+iX)\'] 6{\-\'), 

^ — ' A-/—* 00 

n=0 

(81) 

where we have used the well-known representation of the (5-function 

w^ sinn(A — A') 

5{\-\')^ lim I > . 82 

n^oo TT (A — A') 
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A. 3 Completeness of the set {Pn} using an integral repre- 
sentation 



Using the integral representation of the 2F1 function in the first line of ( |78D and 
the fact that P„ is real, Pn{\,b)* = P„(A, 6), if A real, we can write, inverting 
summation and integration. 



X: PniX, b) P„(A', br = ^T^^^ /o /o dt [s{l - s)t{l - t)r^ (^) (^) ' 



The sum under the integral signs can be evaluated: 

00 

^[(26)„/n!] [(l-2s)(l-2t)]" = [l-(l-2s)(l-2t)]-2b ^ [2{s + t-2st)]- 



2b 



Making the variable substitutions s = e^/(l + e^), t — 6^/(1 + e^) we get for 
the double integral in (p3|), which we denote by / 

poo poo r (x+y)/2-\'^'' 

7 = 2-2''/ / _ e''-^''-^yUxdy. (84) 

J —00 J —00 . ^ ~^ . 

We can "separate" the double integral by the substitution x = u + v, y = 
—u + V to get 

(coshu)^" r(26) ^ 

(85) 

where we have used the integral ([5], p. 11, (26)) 



00 gi2A« 



du = 2 ——rr;du=2^''-'\r{b+iX)\''/r{2b), (86) 



(coshu)^'' Jq (coshu)^'' 

and the representation 

roo 

e'(^-^'> dv ^2TrS{X- X'), 



and furthermore used the ^-function to put A' = A in the first integral. 

Putting (|8^), (|^, and (|5|) together we have obtained the completeness 
relation 

C30 

^ P„(A, b) P„(A', fe) = (5(A - A')/Pfc(A). (87) 

ra=0 
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A. 4 Orthonormality of the set {Pn} using an integral rep- 
resentation 

We want to show that/^ Pm(A, &) P,i(A, 6)pb(A) — 6mn- For this purpose 
it is convenient to use the following contour integral representation of a 2F1 
function with the first parameter equal to a non-positive integer -n, i.e. a 
polynomial of degree n.. The validity is established by expanding the factor 

{1-tz)-''. 

2Fii-n,b,c;z) = \! i f-^'-^l - ty+'^-^l - tz)-'' dt, 



(c)„ 27ri 

where the contour only circles the singularity t = 0. Using ( p8| ) in the orthonor- 
mality integral and inverting integration orders we get 



F„(A, by P„(A, 6) pt{X) dX = v/m!n!/(26),„(2&)„i"+"(-l)"22''-V7rr(26) x 
j jdsdtl^p |r(6 + zA)|2e'^P"(i-2^)"'"(i-2t)]^^| X 



[(l-2.)(l-2t)]'' 

The A-integral in this formula is the Fourier inverse of the integral in ([86|): thus 

/oc 
\T{b + zA)pe'^" d\ = TT2-^^+^T{2b)/{coshu/2f\ (90) 
-00 

Withu = ln[(l-2s)/(l-2i)] weget (coshw/2)2 = {I- s-tf /{l-2s){l-2t), 
and making the coordinate transformations a — s/(l — s), t — t/{l — t) we 
arrive at 

Imn = Vm!n!/(26)„(26)„i"+"(-l)"^^ j j da dr a-"'~^T'''-\\-aT)-^\ 

(91) 

Expanding (1 — err) = X]fe°=o[(2^)fc/^'](''^''')'^ ™^ ^^^^ integral and perform- 
ing the contour integrations we arrive at the desired result 

P„,(A,6)P„(A,6)pb(A) d\ = 5 

ran- 



A. 5 Some comments on the first derivation of complete- 
ness and its relation to Jacobi matrices and the clas- 
sical Hamburger moment problem 

It follows from the discussion in Section 2 , comparing (12), (17), (18), and (31), 
that the completeness relations for the PoUaczek polynomials are (for 6=1/4 
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and 3/4) generalized orthonormality relations for the generalized eigenfunctions 
of the generator Ji = [a^"^ +a^) / A. Correspondingly the set of recursion relations 
761) is the eigenvalue equation Af = Xf for the infinite Jacobi matrix 



/ 



A 



Co 





V 



Co 



Cl 







Cl 



C2 



C2 





Pi 

V •• J 



we have all Cm > 



(92) 



Our case is a special kind of Jacobi matrix in the sense that the diago- 
nal matrix elements of A are all zero; for a general Jacobi matrix they can 
take any (real) values. The connection between the Jacobi matrix and the 
Hamburger moment problem is as follows: The Hamburger moment problem is 
to find a (positive) measure r on the real line corresponding to the moments 
Sm = /_oo u™'dT(u), TO = 0, 1, 2, ... A ucccssary and sufficient condition for the 
existence of a solution is that the sequence {sm} is positive, i.e. that the Hankel 
forms X]r/c=o ^ non-zero vectors {xq,xi, ..} . This condition is 

equivalent to the requirement that a sequence of determinants of the Hankel 
forms should be positive. Certain combinations of the determinants determine 
the coefficients in a Jacobi matrix; expressed in properties of the associated 
Jacobi matrix the condition for the solvability of the corresponding moment 
problem is that all > 0. 

The Jacobi matrix belongs to type D (limit point case; the corresponding 
moment problem is determined, i.e. there is a unique solution for the mea- 
sure), provided ^1/cm — oo (sufficient but not necessary condition!). Then 
the symmetric operator A, defined e.g. on finite sequences /, has a closure 
which is self-adjoint. It is easily seen that the condition on {c™} is satisfied 
in our case. The unique measure is then supported by the spectrum of the 
self-adjoint operator, and is identical to our weight function pb The moments 
Sm — u"^dT{u) can then evidently be obtained from the Fourier transform 
of the weight function (BG) as 



X^"'pb{X) dX = (-l)™[d^"Vdx^"(cosha;/2) 



-261 



(93) 



The odd moments all vanish, corresponding to the fact that the diagonal Jacobi 
matrix elements are zero. 

We refer to ||^ for a more extensive discussion of the connection with Ja- 
cobi matrices and the Hamburger moment problem, in particular regarding the 
case of "higher power squeezing generators" = a'' + a^*^, fc > 3. In each of 
the subspaces = linear span of {\km + k); m = 0, 1, ..}, k = 0, 1, .., A: — 1 
which reduce the eigenvalue equation for Ak the Jacobi matrix derived from 
the recursion relation corresponds to an undetermined moment problem. There 
is an infinity of solutions, and correspondingly the generator Ak is symmetric 
with deficiency indices (1,1) in each subspace. There is in such a subspace a 
one-parameter family of self-adjoint extensions. 
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